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We propose a general method for tailoring the current distribution of relativistic electron bunches. 
The technique relies on a recently proposed method to exchange the longitudinal phase space emit- 
tance with one of the transverse emittances. The method consists of transversely shaping the bunch 
and then converting its transverse profile into a current profile via a transverse-to-longitudinal 
phase-space-exchange beamline. We show that it is possible to tailor the current profile to follow, 
in principle, any desired distributions. We demonstrate, via computer simulations, the application 
of the method to generate trains of microbunches with tunable spacing and linearly-ramped current 
profiles. We also briefly explore potential applications of the technique. 

PACS numbers: 29.27.-a, 41.85.-p, 41.75.Fr 
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I. INTRODUCTION 



Modern applications of particle accelerators generally 
rely on phase space manipulations within one or two 
degrees of freedom. These manipulations are necessary 
to tailor the phase space distributions to match the re- 
quirements of the front-end applications. Usually, phase 
space matching is done in a root-mean-square (rms) 
sense. However, there has been an increasing demand 
for more precise phase space control. In particular, elec- 
tron bunches with tailored temporal distributions are of- 
ten desired. For instance, in accelerator-driven short- 
wavelength light sources based on the free-electron laser 
(FEL) principle and using the so-called high gain har- 
monic generation (HGHG), a uniform current density is 
needed [l| . On the other hand the production of trains 
of microbunches with a given bunch-to-bunch separa- 
tion is foreseen to pave the road toward compact light 
source operating in the super-radiant regime at wave- 
length comparable or larger than the typical microbunch 
length [2|. Finally, novel accelerator concepts based on 
beam-driven acceleration mechanisms^ e.g. plasma or di- 
electric Wakefield- based acceleration iSnOj], would greatly 
benefit from linearly ramped or trapezoidal current pro- 
files [3, Q to significantly increase the transformer ra- 
tio - the energy gain of the accelerated bunch over the 
energy loss of the driving bunch [3[. Furthermore the 
performance of the aforementioned acceleration schemes, 
could be significantly improved if driven by a train of 
microbunches to resonantly drive wakefields [9| or train 
of linearly ramped microbunches jlOl . Ill| to improve the 
transformer ratio. 

To address these needs, several techniques aimed at tai- 
loring the current profile of electron bunches have been 
explored. Linearly-ramped current profiles can be pro- 
duced by imparting nonlinear distortions in the longitu- 
dinal phase space |12h14| . A successful experimental im- 



plementation of a nonlinear optical lattice was reported 
in Ref. [15|. The generation of a train of microbunches 
by shaping the temporal profile associated to a photo- 
cathode drive laser in a photoemission electron source 
were explored via numerical simulations |16l4l8| . An al- 
ternative method using an interceptive mask located in 
a dispersive section, first proposed in Reference [l9| , was 
experimentally demonstrated [20]. Tailoring of the cur- 
rent profile to enhance the operation of a FEL operating 
in the HGHG regime via an ad hoc temporal shaping of 
the drive-laser was discussed in Reference [21| . 

Each of these methods has limitations. The techniques 
based on shaping of the photocathode drive laser distri- 
bution are prone to the space-charge effects which are 
prominent at low energies and tend to smear the im- 
pressed shaping. Although it has been pointed out that 
the density modulation is indeed converted into an en- 
ergy modulation and can be later recovered by a proper 
longitudinal phase space manipulation, the final bunch- 
ing factor is significantly lower ^^. Similarly the tech- 
niques based on nonlinear transformations generally have 
a limited tunability. 

In this paper we propose an alternative method to 
tailor the current distribution of relativistic electron 
bunches extending on our earlier work [2j, [2J|- The 
technique relies on the transverse-to-longitudinal phase 
space exchange technique. After discussing the theoreti- 
cal background for our proposal in Section [Hi we demon- 
strate the generation of various current profiles via com- 
puter modeling in Section Hill We finally discuss a pos- 
sible implementation and its applications in Section [IVI 



II. THEORETICAL BACKGROUND 

The proposed method for producing beams with ar- 
bitrary current profile is illustrated in Figure [T] A 
transversely-shaped electron bunch with the transverse 



density distribution ^±{x,y) is manipulated through 
a transverse-to-longitudinal emittance exchange (EEX) 
beamline to convert the corresponding horizontal profile 
P{x) — J dy^±{x, y) into the longitudinal (temporal) co- 
ordinate; see Fig. [TJ 
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FIG. 1: (Color) Overview of the proposed technique to pro- 
duce relativistic electron bunch with arbitrary current profile. 

The backbone of our proposed scheme is the 
transverse-to-longitudinal phase space exchange. The 
emittance exchange was first proposed in the context of 
B-factory [25] as a way to achieve very small /3* values 
at the interaction point. The scheme was later explored 
as a possible alternative for mitigating microbunching 
instability in bright electron beams [26] or for improv- 
ing the performance of single-pass FELs [27[ . There are 
several solutions capable of performing this phase space 
exchange [2^ [23| ■ The simplest solution devised to date 
consists of a horizontally-deflecting cavity, operating on 
the TMiio mode, flanked by two identical horizontally- 
dispersive sections henceforth referred to as "dogleg" [30| ; 
see Fig. [TJ An experiment conducted at the Fermilab's 
AO photoinjector facility recently demonstrated the ex- 
change of transverse and longitudinal emittances [3i| . 

We first elaborate on the emittance exchange princi- 
ple and analyze the single-particle dynamics in the four 
dimensional phase space (x,x',z,^), an initial electron 
with coordinate Xp = (xq, Xq, zq, ^o) is mapped to its fi- 
nal coordinate X = (z, a;', z, 5) accordingly to X = fi Xp 
wherein R is the 4x4 transport matrix; the ~ stands 
for the transposition operator. In our notation under- 
lined vectors belong to the four-dimensional trace space 
(x, x' , z, 5). The matrix of a dogleg composed of two rect- 
angular bends of length Lb and bending angle a sepa- 
rated by a drift space D has the form 
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where the matrix elements are related to the machine 
parameters by the following: L = l£°Jq+Il^ ^ = 



and r/ 



D+2c^Li-cl(2Lb + D) 



with 



the short-handed notation {c^s)a = (cos, sin)a [32]. The 
transfer matrix associated to a deflecting cavity of length 
Lc and deflecting strength n is [2a, [33] : 
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A complete emittance exchange requires k 
and the total matrix becomes: 



M = MdMcMd 



A B 
C D 



(2) 
(3) 



where the elements of matrix A are all zero except A12 

{Lc+AL)/4: (4^2. - 4772 +^L,)/4 

has the same form as B but its diagonal elements are re- 
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versed, and D = j^ 



1 ? 



. For a thin cavity {Lc = 0), 



the matrices A and D vanish and the exchange is per- 
fect [30] • In this latter case we have a direct mapping of 
the transverse to longitudinal coordinates and vice versa. 
In particular we have: 



(4) 



In general, when including thick lens effects, the final 
beam matrix is E = ( XX ) = MSpM is not block anti- 
diagonal and the final geometric emittances are given to 
first order by [2^ 
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where qq = Sxfi/£z.o is the initial transverse-to- 
longitudinal emittance ratio and A^ is given by [32| 



A^ = 



64.rj^Px.oPz.o 



(6) 



Here a^p and /3j are the Cour ant- Snyder (C-S) pa- 
rameters associated to the horizontal {i = x) and 
longitudinal (i = z) degrees of freedom. The quantity 
A^ can be minimized by a proper choice of either 
longitudinal or transverse C-S parameters. Here we 
opt for minimizing A with respect to the longitudinal 
C-S parameters. Since the bunch length (i.e. Pzfl) is 
fixed by the upstream beamline and electron source 
settings for most of the cases discussed in this paper, A 
is minimized with respect to az q = — ° ° , the opti- 
^z,o = f corresponding to an incoming 
'^'^'1 - -i that 



mum value is az 

longitudinal phase space chirp C = ||] 



produces a minimum bunch length at the cavity location. 

In order to illustrate how this phase space exchange 
mechanism could be taken advantage of to tailor the lon- 
gitudinal distribution of an electron bunch, we consider 
a beam with an initial six-dimensional Gaussian trace 
space distribution 



*o(Xo) 



no 



1. 



/o ^3=- . — exp[--XoE/Xo] (7) 



where Xq = {xo,x'o,yo,yo,zo,So), and S = (XqXo) 
is the six-dimensional covariance matrix associated to 
the beam. We further assume the initial covariance 
matrix to be uncoupled, i.e. to have the form Eq = 
diag{T,a;fi,^v,o,^z,o) with 



Ei,0 = £i,0 



A,0 —Oii.Q 



(8) 



where the index i = x,y, z stands for the considered de- 
gree of freedom and 7i^o = (1 + oifoj/l^ifi- 

We suppose the initial distribution function can 
be transversely shaped with a "transmission" function 
T(xo,yQ). The function T(a:o,2/o) can model a mask [if 
so its value is either 1 (within the mask aperture) or 
(elsewhere)] or be used to account for a transverse density 
that deviates from the Gaussian distribution. Therefore 
the four-dimensional trace space distribution upstream 
of the EEX beamline is 



<&o+(Xo) 



'&o~(Xo;0) 



X exp 



4-00 



dyo 



vl 
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Downstream of the EEX beamline the final distribution 
is of the form 



$(Xo;s) - $(i?-iXo;0)T(i?-iXo 
and the longitudinal trace space is given by 
no 



(10) 
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2ex,a 

XTiNi5z + NieS), 

are the elements of the matrix N 



(11) 

The latter equation describes a two-dimensional Gaus- 
sian distribution in {z,S) modulated by the function r. 
Integration over 6 provides the longitudinal profile down- 
stream of the exchanger: 



P{z) = 



no 



2TTe 



x,a 



■ exp 



'^l^z^xfi 
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with Q{z) — J_ exp 



P. 



13. 



r{Ni5Z + 



NieS)dS. If the mask has no y-dependence [y^{x) -^ +oo 
and y^ (x) — >■ — cxd ], r is a function of x only. 

We first consider the case of t(xo) = J2i=i ^i^ ~~ ^o) 
where 6{...) is the Dirac function. This can be achieved, 
e.g., by intercepting the incoming beam with a series 
of vertical slits with center located at Xq. In such a 
case the longitudinal phase space is described by a two- 
dimensional Gaussian distribution with C-S parameters 
(azjfSz), emittance e^o, which is non- vanishing only at 
locations described by the family of lines with equation 
6{z) = {xq + Ni5z)/Ni(i in the longitudinal phase space; 
see Fig. [2] The associated projection along the longitu- 
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FIG. 2: Rendition of a stratified longitudinal phase space 
generated downstream of the EEX beamline with an incoming 
transversely-modulated beam. 

dinal coordinate z is 
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representing a series of Gaussian microbunches each 
of length (7z,7n = \^£x,oNiq/I3z. The microbunches 
are longitudinally separated when 4:az^m ^ Xq/[Ni^ — 
<^z/ PzNie]. Even when the P{z) projection is not 
strongly modulated, the local correlations in the longitu- 
dinal trace space can be taken advantage of to enhance 
the modulation. This can be done by using a disper- 
sive section downstream of the EEX beamline with its 
longitudinal dispersion ^' matched to the chirp of the 
beamlets, i.e. £,' = —Ni^/Niq. A simple estimate of ^' in 
the thin-lens approximation gives the required value to 



be 



i'-i-\><^^ 



(14) 



In our convention [head of the bunch is at positive time 
(or z-coordinate)] the required value is of opposite sign 
of the value provided by, e.g., a standard four-bend 
magnetic chicane. Such a longitudinal dispersion can 
be generated by a four-bend chicane with quadrupoles 
inserted between the two first and two last dipoles |34| . 

We now consider the more general case of a mask with 
upper (y > 0) and lower {y < 0) boundaries respec- 
tively respectively described by the functions yQ{xo) and 
?/(^(xo). In such case we have 



t(Xo) = T{xa) = j 



X exp 






dyo 



T{xQ,yo). 



(15) 



Taking j/q — —y^ = (^(xo), where the function ^{x) 
describes the mask shape, and T{xq) = 1 for ya{xa) G 
[yo{xo),y^{xa)] yields 



'[xo) — 2 erf 
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(16) 



The general form of the longitudinal bunch profile 
downstream of the EEX bcamline is 
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The latter equation provides some insight on the limi- 
tation of the discussed pulse shaping technique. For in- 
stance, the convolution-like integral indicates that cur- 
rent profiles with spectral components beyond the spa- 
tial frequency fx ~ l/[2Try^/3z£x,o] will be filtered out. 
Therefore the generation of pulses with sharp features 
would require the final longitudinal emittance and beta- 
tron function to be small. 



III. VALIDATION OF THE METHOD VIA 
SINGLE-PARTICLE DYNAMICS SIMULATIONS 

We first validate the technique proposed in the previ- 
ous Section with the help of single-particle dynamics sim- 
ulations. We consider an incoming beam with Gaussian 
distribution in the three degrees of freedom. The beam 
is intercepted by a mask with transverse transmission 
T(a:, y) and then passes through the emittance exchanger. 
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FIG. 3; Evolution of the transverse horizontal (black), ver- 
tical (cyan) and longitudinal (red) emittances along the 
EEX beamline with parameters given in Table. [H The 
green symbols indicate the axial locations of the four dipoles 
(shown as rectangles) and deflecting cavity (vertical line at 
s ~ 3 m). The initial normalized emittance partition is 
(eS,o,ea,o,£"o) = (1,1,10) ^m. 



The beamline geometric parameters are gathered in Ta- 
bleland correspond to the proof-of-principle experiment 
being considered at the Argonne Wakefield Accelerator 
(AWA) facility [35|. The transversely-deflecting cavity 
consists of a 1/2-1-1/2 cell cavity optimized to minimize 
steering effects [36| . The numerical beam dynamics simu- 



TABLE I: Geometric parameters and properties associated to 
the EEX beamline considered in this paper; see also Fig. [T] 



parameter 



symbol value unit 



Bending angle 

Bend efi'ective projected length 

Projected distance between bends 

Dogleg final dispersion 

Dogleg longitudinal dispersion 



a 


20 


deg 


Lb 


0.308 


m 


D 


0.97 


m 


n 


0.51 


m 


i 


0.16 


m 



lations presented in this section were performed with the 
particle tracking code elegant ^37] and do not include 
any collective effects. The influence of collective effects is 
investigated in the next section when considering a realis- 
tic beamline configuration along with exploring practical 
applications of the pulse shaping technique. 

In ELEGANT the deflecting cavity was modeled by its 
thick-lens transfer matrix which was found to be in ex- 
cellent agreement with the transfer matrix numerically 
devised by tracking macroparticles in the 3D electromag- 
netic field map computed with the microwave studio® 
eigensolver |38|; see Appendix A. An example of emit- 
tance evolution along the emittance exchanger with pa- 
rameters tabulated in Table. HI is shown in Fig. |31 



Limiting effects in an emittance-exchange 
beamline 



downstream of the mask is then used in elegant to 
model the emittance exchanger beamhne; see Fig. [S] 



Before proceeding with vahdating the pulse shaping 
method described in the previous section, we explore the 
domain of perfect phase space exchange. Several sources 
can prevent a perfect emittance exchange. One of them, 
described in the previous section, is due to coupling in- 
troduced by the thick lens model of the deflecting mode 
cavity. Although minimizing this emittance dilution re- 
quires the chirp to be tuned to minimize the bunch length 
at the cavity location, the emittance is still diluted and 
the dilution term can be further minimized by select- 
ing proper incoming horizontal C-S parameters as illus- 
trated in Fig. m In the latter Figure we introduced the 
normalized emittances defined for relativistic beams as 
e"o = 7ei,o where i ^ x,y,z. 




FIG. 4: Fractional horizontal (left column) and longitudi- 
nal (right column) emittance dilutions at the exit of the ex- 
changer beamline versus the incoming horizontal Courant- 
Snyder beam parameters. These quantities are respectively 
defined as Se^ = ex/£z,o — 1 and Ssz = £z/£x.o ~ 1- The 
rows correspond to dilTerent incoming normalized emittance 
partitions from top to bottom (£",oiS"o) = (IO,l), (1,1) and 
(1,10) [all numbers are in micrometers]. The longitudinal 
phase space parameters used in these calculations are a bunch 
length (Tz,o — 400 ^m, and a chirp C — % m~^. 
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FIG. 5: Example of conversion of a transverse modulation (a) 
into a modulation on the longitudinal phase space (d). The 
plots (a) and (c) show the initial horizontal and longitudinal 
phase space while the plots (b) and (d) show the correspond- 
ing final phase spaces downstream of the EEX beamline. The 
initial emittance partition is (e" q, e"o) = (l,10) [ra\. 



A series of simulations were performed for three cases 
of incoming transverse-to-longitudinal emittance ratio 
Po- The horizontal C-S parameters were chosen for an 
optimum emittance exchange in accordance with Fig. |4] 
and the beam's Lorentz is 7 = 30. The obtained cur- 
rent profiles were characterized by computing the modu- 
lation wavelength Am , and the associated bunching factor 
hra = b{km) whcrc 



5(fc)- 



+ 00 



P{z)e'''^dz 



and km — 27r/Am- The contrast ratio 



X = 



P -P 

P + P' 



(18) 



(19) 



B. Generation of stratified longitudinal phase 
spaces and train of microbunches 

We first demonstrate the generation of a train of mi- 
crobunches using an interceptive mask located upstream 
of emittance exchanger beamline. For these simulations 
the incoming Gaussian beam is intercepted by a 3 mm 
thick tungsten mask. The slits width is taken to be 
100 fim. and the edge-to-edge separation between the 
slits, Sw, is variable. The simulations of the beam 
interaction with the mask, including scattering effects, 
is executed with shower j39j and the beam distribution 



where P = max[P(z)] and P = min[P(z)], was also cal- 
culated. The resulting values for Am, &m and the cor- 
responding contrast ratio Xm are presented in Figure El 
We also include the parameters computed after removal 
of the local correlations such to result in upright lon- 
gitudinal trace space for the beamlets: the modulation 
wavelength Xm, associated bunching factor 6jv/ = h'{kM) 
(where kM = Stt/Aa/ and b' is the bunching factor com- 
puted after the longitudinal phase space manipulation), 
and the contrast ratio xm ■ When the final horizontal and 
longitudinal emittances are equal [(e" q, s" o)=(1'10) M^i] 
or when the final longitudinal emittance is smaller than 



the horizontal one [(e" g, e" g)=(l,10) ^m], a clear bunch- 
ing is observed (for most of the cases b^ > 0.1 and 
Xm ^ 0.5). For the case when the final longitudinal 
emittance is larger [(e" Q,e"g)=(10,l) /.tm], the bunch- 
ing factor is small, i.e. bm ^ 10~^ but the longitudinal 
trace space is clearly stratified and the removal of the 
local correlations leads to strong bunching: the bunch- 
ing factor becomes Bm ^ 0.5 and the contrast ratios for 
all cases are close to unity. The modulation wavelength 
for the considered cases was 0.1 < Xm < 0.2 mm and 
0.035Am < 0.080 mm. 



parameters (a^.o and /?s.o) upstream of the exchanger 
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FIG. 6; Bunching factor (a), contrast ratio (b) and modula- 
tion wavelength (c) obtained for a longitudinally modulated 
beam produced via interception with thick mask upstream of 
the EEX beamline. The slits width is w = 100 /im and the 
slits separation Sw is variable. The parameters with subscript 
m (and shown as blue traces) are computed using the longi- 
tudinal beam distribution downstream of the EEX beamline 
while the parameters with subscript m (and shown as red 
traces) are obtained after the removal of the chirp associated 
to the beamlet (in this latter case, the beamlets are made 
upright) . The parameters are calculated for various initial 
normalized emittance partitions (e",oi£"o) shown in the leg- 
end. The beam Lorentz factor for these simulations is 7 = 30. 



Ps,0 = "^Q^s.o ± 



^) -z/L^?Kl + <o) 



1/2 



, (20) 



where v = {1 — rf)/{L^). For f3sQ to be real-positive the 
condition |as_o| > |1 + v\l \fv needs to be fulfilled. Equa- 
tion [20] represents a hyperbola in the (as.Oj/3s,o) space. 
The horizontal C-S parameters associated to a beamlet 
generated by intercepting a beam with initial parame- 
ters (a_,/3_) are as,- = 0, /3s__ = /3_ and the beam- 
let's emittance is £«,_ = w^/(12/3_) where w is the slit 
width. Therefore a set of quadrupoles would be needed 
between the mask and the EEX beamline to match the 
initial beamlet C-S parameters (q;_,/3_) to the required 
values (as,o,/3s.o) upstream of the EEX beamline. Do- 
ing so will also modify the C-S parameters of the entire 
beam and might result in a significant beam emittance 
dilution as discussed in Section IIII Al For the geome- 
try considered in this paper (see Tab. H]) and under the 
single-particle dynamics it is possible to tune the C-S 
parameters to get the beamlets' longitudinal trace space 
upright downstream of the EEX beamline without sig- 
nificantly diluting the emittances; see Fig. [T] In general, 
it is always possible to find a beamline configuration up- 
stream of the EEX beamline to simultaneously provide 
upright trace space for the beamlets while matching the 
horizontal C-S parameters of the entire beam to values 
minimizing emittance dilution. This can be done by con- 
structing a beamline that provides the following incoming 
beam parameter upstream of the slits: 



a_ 







^/K^s 






(21) 



Therefore inserting the mask between two matching sec- 
tions, each minimally composed of four quadrupoles, en- 
ables to control the beamlets longitudinal trace space cor- 
relation downstream of the EEX beamline while minimiz- 
ing beam emittance dilution during transport through 
the EEX beamline. 



In Section II, we briefly discussed the possible use of 
a properly designed dispersive section to remove the cor- 
relation and produce a train of microbunches. In fact, 
another scheme is to choose the proper initial horizon- 
tal C-S parameters upstream of the EEX beamline such 
that the final longitudinal trace spaces associated to each 
beamlet are uncorrelated downstream of the EEX beam- 
line. For the sake of simplicity we consider a simple model 
for the EEX beamline where the deflecting cavity is de- 
scribed by its thin-lens approximation; see Eq. |4l Im- 
posing the flnal longitudinal trace space associated to, 
e.g., the central beamlet to be uncorrelated [(zS) — 0) 
yields a condition between the beamlets horizontal C-S 



C. Arbitrary current profiles 

The generation of arbitrarily-shaped current profiles 
could be achieved using a mask technique as described in 
the previous section. As an illustrative example, we first 
consider the academic case of a transverse density follow- 
ing a uniform triangular distribution; see Fig. [5] (a). Such 
a distribution results in a linearly-ramped current profile 
downstream of the EEX beamline upon proper choice of 
C-S parameter; see Fig. [8](d,e). The latter figure clearly 
demonstrates the inability of the pulse shaping technique 
to resolve the initial high-frequency features (in this case 
the hard edge associated to the triangular profile). As 




FIG. 7: Locus of C-S parameters that provide an upright 
longitudinal trace space for the beamlets (white solid line). 
The false color map represents the emittance dilution 5tz as 
a function of initial C-S parameters and the dashed white line 
indicate the limit where emittance dilution is 5%. The initial 
normalized emittance partition is (£",oi£z,o) = (lOi 1) A*™- 
The solid white line is computed numerically and corresponds 
to the Eq. [20] 



discussed in Section|TT]this limitation stems from, the final 
longitudinal emittance values. In the present case, for the 
same final longitudinal betatron function, the beam with 
larger final longitudinal emittance has its temporal pro- 
file significantly smeared [see Fig. [S^d)] compared to the 
beam with lower longitudinal emittance [see Fig. [8l^e)]. 

For realistic beams , a mask with suitable shape needs 
to be used. Taking the case of a initial Gaussian hori- 
zontal profile and imposing the horizontal profile down- 
stream of the mask to be a linear function of z [P{x) oc x\ 
results in the mask shape 



v{x) = v2crj,erfinv 



' 51 xj exp 



2(Ti 



, (22) 



where g^ and gi are constants that control the offset and 
steepness of the linear ramp. The mask profile along 
with the initial and final horizontal profiles are depicted 
in Fig. [9] For this particular case of an initial Gaussian 
distribution, only ~ 19% of the incoming charge is trans- 
mitted through the mask thus making the masking tech- 
nique highly inefficient. The case of an initial uniform 
circular distribution is more favorable but still results in 
excessive charge loss. 

In principle any current profile can be synthesized with 
the proposed technique. The main disadvantage of using 
a mask being the substantial charge loss. An alterna- 
tive method is to shape the beam using the technique 
discussed later in Section HVBI 
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FIG. 8: Example of the generation of linearly-ramped cur- 
rent profile from an initial uniform triangular distribution in 
the (x, y) space (a) shown with corresponding final transverse 
distribution (b) and longitudinal phase space (c) downstream 
of the EEX beamline. The initial normalized emittance par- 
tition is (e"_o, e" o) = (10, 1) Mm and {a^,o, Px,o, ) = (2, 12 m). 
The current profiles downstream of the EEX beamline are 
shown in plots (d) and (e) for respectively {e"fi, e" o) = (10, 1) 
and (1,10) ^m. The incoming horizontal C-S parameters are 
a betatron function of 12 m and the a-function was varied 
by increment of 1 from -4 to +2 [corresponding to the seven 
traces shown in plots (d) and (e)]. The initial beam Lorentz 
factor for these simulations is 7 = 30. The head of the bunch 
corresponds to f < 0. 



IV. POSSIBLE APPLICATIONS IN 
LOW-ENERGY ACCELERATORS 

In this Section we demonstrate, via start-to-end sim- 
ulations, possible use of the current shaping technique 
described in this Paper. We use as a simulation platform 
the AWA facility shown in Fig. [TU] The accelerator in- 
corporates a photoemission source consisting of a 1-1-1/2 
cell rf cavity operating at / = 1.3 GHz, henceforth 
referred to as rf gun. An ultraviolet (uv) laser beam 
impinges a magnesium photocathode located on the back 
plate of the rf gun half cell. The thereby photoemitted 
electron bunch exits from the rf gun with a maximum 
kinetic energy of approximately 8 MeV and is injected 
into a 10-cell standing wave normal conducting booster 
cavity operating on the 7r/2 mode at / = 1.3 GHz. 
The final beam energy can be up to ~ 16 MeV. The rf 
gun is surrounded by three solenoidal magnetic lenses 
independently powered referred to as LI, L2 and L3; see 
Fig.Iini 

The simulations of the EEX beamline presented in this 
Section were performed with the particle-in-cell (PIC) 
program impact-t which includes space charge effects by 
solving Poisson's equation in the bunch's rest frame us- 
ing a three-dimensional quasi-static algorithm [40l|. The 
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FIG. 9: Example of mask shape needed to generate a linearly- 
ramped horizontal profile by interception of beam with trans- 
verse radial Gaussian density (a). The mask profile used is 
described by Eg. 1221 and only macroparticles with \y\ < vix) 
are kept resulting in the transverse charge density shown in 
(b). The initial Gaussian horizontal profile (c) is transformed 
into a linear-ramped horizontal profile (d). For this exam- 
ple the initial transverse beam size are a^ ~ 1.4 mm and 
ay — 2.5 mm and the parameters in Eq. [22]are go = 100 and 
ffi = 3 X 10''. 



(BFF). Considering a series of iVf, identical microbunches 
with normalized distribution A(t) we have S{t) = 
N-^ J2n=i Ht + nT) (where T = A^/c is the period and 
Am the density modulation wavelength defined earlier) 
giving |S'(a;)p — E\A{uj)\^. The intra-bunch coherence 
factor S EE 7Vj"2 sin'^ {coNbT/ 2) /[si-a^icoT/ 2)] describes the 
enhancement of radiation emission at resonance, i.e. for 
frequencies w = 27rpc/A„i (where p is an integer). 

In order to demonstrate the generation of a train of mi- 
crobunch or modulated bunches using realistic distribu- 
tions produced by a photoinjector, we modeled the AWA 
photoinjector from the photocathode up to downstream 
of the booster cavity with ASTRA (up to z = 2.79 m in 
Fig. [ini), a PIC program that incorporates a cyhndrical- 
symmetric space charge algorithm 4q|. The generated 
1-nC bunch was coUimated through a set of vertical 
slits thereby producing beamlets and its final horizon- 
tal C-S parameters upstream of the EEX beamline were 
matched to specific values using a set of four quadrupoles 
shown as Ql, Q2, Q3, and Q4 in Fig. ^U\ The track- 
ing from z — 2.79 m up to z ~ 6.4 m was performed 
with IMPACT-T and includes space charge effects. The iV- 
macroparticles ensemble downstream of the EEX beam- 
line was represented as a Klimontovich distribution 



deflecting mode cavity was simulated with microwave- 
STUDIO® and the generated 3-D electromagnetic field 
maps were imported in impact-t; see Appendix A. The 
simulated magnetic field of the horizontally-bending rect- 
angular dipoles was used to accurately model the dipoles 
in iMPACT-T by representing the fringe field regions with 
eight-order Enge functions [41 1. 



A. Terahertz radiation source 

Terahertz (THz) radiation occupies a very large por- 
tion of the electromagnetic spectrum and has generated 
much recent interest due to its ability to penetrate deep 
into many organic materials without the damage associ- 
ated with ionizing radiation such as x-rays. THz radia- 
tion is easily produced in laser-based method via optical 
rectification |42l444| . The radiation spectrum is generally 
broadband and the associated energy is of the order of 
nJ per pulse. Accelerator-based THz radiation sources 
present advantages in term of emitted power, frequency 
tunability and spectral bandwidth. The spectral fluence 
radiated by a bunch of N^ electrons taken to follow a 
line-charge distribution is j45i] 



dP 

dcddQ 



We 



dP 

dujdfl' 



]i N, + N,{N,~1)F{lj) (23) 



where [ jZdn ] i ^^ ^^^ single electron power density as- 
sociated to the considered electromagnetic process and 



1 ^ 



Zj) 



(24) 



along the longitudinal coordinate. The associated BFF 
is calculated as the summation 



F{uj) 



N 

E 



sm{ujZj/c 



N 



^cos{ujzj/c) 



i=i 



(25) 



Figure [TT] depicts examples of modulated bunches (the 
charge is 300 pC downstream of the slits). The BFF 
associated are enhanced at the fundamental bunching 
and harmonic frequencies which can be varied by either 
changing the quadrupoles settings upstream of the EEX 
beamline or using different slits separations; see FigfT^ 
In this particular example the frequency can in princi- 
ple be continuously varied from ~ 0.6 to ^ 2.2 THz (if 
we consider the fundamental bunching frequency only; 
higher frequencies can be obtained by proper filtering to 
select harmonics of the bunching frequency) . We should 
note that the calculated BFF is representative of the ra- 
diation spectrum as long as the line charge distribution 
assumption is fulfilled, e.g. provided the longitudinal and 
transverse a± rms bunch sizes satisfy az ^ cr±/"f if the 
considered radiation process is transition radiation (47| . 
The generated train of microbunches can also be used to 
produce super-radiant undulator radiation [2|. 
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FIG. 10: Overview of the photoinjector configuration used for reahstic simulations of bunch shaping. The legend is "Di" dipole 
magnets, "Qi" quadrupole magnets, "Li" solenoidal lenses. The geometry of the EEX beamline is detailed in Table HI 
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FIG. 11: Example of conversion of modulated longitudinal 
trace space downstream of the EEX beamhne. The four plots 
corresponds to the four cases of initial parameters (SwjWjT) 
of (100, 300, 0.92) [plot (a)], (100, 300, 1.25) [plot (b)], 
(50, 150, 0.75) [plot (c)] and (25, 75, 0.75) [plot (d)] where 
Sw and w are respectively the slit edge-to-edge separation 
and width in micrometers, and J^ is a. scaling factor for 
the quadrupoles Ql, Q2, Q3, and Q4 (when J^ = 1 the 
quadrupoles are tuned to match the horizontal C-S param- 
eter to {ax,o,Px,o) = (2,12 m) upstream of the exchanger). 



B. Ramped drive beams and bunch trains for 
beam-driven wakefield accelerators 




10 10 

f (THz) 

FIG. 12: Bunch form factor associated to train of mi- 
crobunches obtained downstream of the EEX beamline. The 
four traces correspond to four settings/configuration indi- 
cated in the legend as "Sw — w — J-" where 5w and w are 
respectively the slit edge-to-edge separation and width in mi- 
crometers, and J^ is a scaling factor for the quadrupoles Ql, 
Q2, Q3, and Q4 (when J^ — 1 the quadrupoles are tuned to 
match the horizontal C-S parameter to {ax,o,l3x,o) ~ (2, 12 m) 
upstream of the exchanger). 



the photocathode. Several types of distribution were con- 
sidered for the initial laser transverse distribution and we 
found that a distribution with a dipole moment in pro- 
vides the desired ramped horizontal profile upstream of 
the EEX beamline. We therefore consider laser trans- 
verse densities on the cathode described by [48| 



Another application of pulse shaping we consider is 
the generation of beams with linearly-ramped current 
profiles. Ramped beams maximize the transformer ratio 
in coUinear beam-driven acceleration schemes. A prac- 
tical implementation using an intercepting mask seems 
highly inefficient as discussed in Section IIII CI Instead 
we explore the generation of such a beam by generating 
a beam with a suitable transverse distribution starting at 






(26) 



where po{r,9) — 1 over a specified domain and else- 
where, Tc is the beam size on the photocathode, ^o is the 
initial rotation angle of the dipole pattern, and A is the 
amplitude of the perturbation. The symmetry axis of the 
laser distribution is rotated such that after propagation 
through the emittance-compensating solenoid the x — y 
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coupling is removed. This is accomplished by choosing 
the rotation angle to be |49| : 



[eB{r ^ 0, z)]/[2Tn-f{z)l3{z)c]dz (27) 



where B[r — 0, z) is the axial component of magnetic 
field experienced by the beam, 7 the beam's Lorentz fac- 
tor and /3 = [1-7-2] 1/2. 

Three types of distributions were considered: (1) an 
initial uniform triangle (po — 1 inside a equilateral trian- 
gle and ^ = 0), (2) a circular distribution with a dipole 
moment (po — 1 for r < re and V0 and ^ = 1), and (3) 
a triangular-shaped distribution with a dipole moment 
{A— 1). For all cases the distribution on the cathode was 
centered w.r.t to its barycenter and the rms transverse 
sizes were forced to be 1 mm. These three types of dis- 
tribution can be implemented using uv masks of variable 
transmission [50| . The laser density on the photocathode 
surface and the resulting 1-nC bunch at z = 2.79 m are 
shown in Fig. [T^ 
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FIG. 13: Initial transverse densities on the photocathode sur- 
face (top row) and corresponding transverse densities (middle 
row) and horizontal profiles (bottom row) at 2; = 2.79 m. The 
three column respectively correspond to the case of a uniform- 
triangle distribution (left column), a circular distribution with 
a dipole moment (middle column) and a triangular-shaped 
distribution with a dipole moment (right column). 

Figure [13] demonstrates that these distributions have a 
ramped horizontal profile upstream of the EEX beamline. 
For the considered accelerator settings, the triangular- 
shaped distribution with a dipole moment yields the best 
linearly-ramped horizontal profile that can be converted 
into a linearly-ramped current profile downstream of the 
EEX beamline as illustrated for a 1-nC bunch in Fig. [TD 
The quadrupoles Ql, Q2, Q3, and Q4 were tuned to ap- 
proximately provide C-S horizontal parameters upstream 
of the EEX beamline close to (2, 12 m) (the matching was 
done in the single-particle dynamics mode and the cor- 
responding quadrupole settings were used in impact-t 
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FIG. 14: Example of longitudinal phase space (left) and as- 
sociated linearly-ramped current profile (b) obtained down- 
stream of the EEX beamline. The beam Lorentz factor is 
7 ~ 33.4 and the rms bunch length is 490 /im. The initial 
distribution used for this calculation is the triangular-shaped 
distribution with a dipole moment (right column in Fig. [13]). 



The final normalized emittance partition (e", e" e") = 
(43.6,5.8,12.1) /im makes such a bunch suitable as a 
drive beam for beam-driven techniques based on asym- 
metric structures, e.g., dielectric slabs [5ll[54|. The final 
current distribution is used to estimate the achievable ax- 
ial Wakefield following Ref . [5i| . We consider a slab struc- 
ture with half gap a = 0.150 mm and dielectric thickness 
h — a = 0.150 mm We take the relative dielectric per- 
mittivity to be e = 9.8 and compute the axial wakefield 
produced by the bunch as 



N 



^-(o = EE"'^(c-o 



(28) 



-1 j=i 



where (^ = z — Vbt {vb ~ c is the bunch's velocity), Wi(C) 
is the Green's function for the i-th mode wakefield de- 
scribed by Eq. 2.6 in Ref. j^]. The indices i and j in 
the latter equation correspond respectively to the sum- 
mations over the number of modes considered (we take 
n = 10) and the number of macroparticles (iV = 2 x 10^) 
representing the bunch (here we model the bunch as a 
Klimontovich distribution; see Eq. [24]) . The computed 
axial electric field is shown and compared with the elec- 
tric field induced by a Gaussian beam with similar rms 
duration in Fig. [T^] . The transformer ratio defined as 
the ratio of the maximum accelerating field to the maxi- 
mum (in absolute value) decelerating field is found to be 
n ~ 20.08/10.22 = 1.96 and 11 ~ 22.05/6.13 = 3.60 for 
respectively the Gaussian and linearly-ramped bunch. 

If a round beam is desired downstream of the EEX 
beamline, the shaping technique would need to be com- 
bined with the fiat beam technique [5J, [5J| . 

Finally, a bunch with linearly-ramped horizontal pro- 
files could also be intercepted with a series of slits to 
generate a train of bunch. This would result in a ramped 
bunch train and has application to resonantly excite 
wakefield in beam-driven acceleration schemes. 
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FIG. 15: Axial electric field (blue traces) behind a 1-nC bunch 
with Gaussian (a) and linearly-ramped (b) current profile (red 
traces). In our convention Ez < corresponds to a deceler- 
ating field. For both case the bunch is centered at ^ = and 
the head corresponds to ^ > 0. The transformer ratio quoted 
in the text is computed as the ratio of Ez{(^ ~ —1.3 mm) over 
|£.(C==;Omm)|. 



SUMMARY 
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APPENDIX A: NUMERICAL VALIDATION OF 

THE ANALYTICAL MATRIX OF A 

DEFLECTING CAVITY 

In this appendix we show that the transfer matrix of 
a deflecting cavity derived in Ref. [23 for an ideal pillbox 
cavity operating on the TMhq accurately describes the 
matrix modeled numerically computed using the electro- 
magnetic field obtained from microwave-STUDIO® sim- 
ulations. 



We presented a general method for generating relativis- 
tic electron beams with arbitrary current profiles. The 
proposed method offers a very generic tool for precisely 
controlling the current distribution (and possibly the 
emittance partitions) to match the requirements imposed 
by the front end application. The technique takes advan- 
tage of the recently proposed transverse-to-longitudinal 
phase space exchange to map the horizontal profile asso- 
ciated to the electron bunch into its current profile. In 
principle any current profile can be produced by prop- 
erly manipulating the bunch's transverse density to pro- 
vide the desired horizontal profile upstream of the EEX 
beamline. We have considered special cases and showed 
how the method could be used to produce a train of 
microbunches with variable spacing and linearly-ramped 
bunch using start-to-end simulations of a realistic con- 
figuration in a ^ 16 MeV photoinjector. The technique 
could also be applied to high-energy beams as discussed, 
for instance, in Ref. [5a|. At high energies, the use of a 
thick mask to intercept the beam might be prohibited and 
instead a selection technique based on emittance spoiling 
using a thin foil might be preferable |56|. 

Note added in proof — Very recently an experiment 
conducted at the Fermilab's AO photoinjector demon- 
strated some of the possibilities discussed in this pa- 
per [57| including the generation of a train of sub- 
picosecond bunches with variable separation [58|. 
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FIG. 16: Geometry of the TMoio-like defiecting mode cavity 
used for the simulation presented in this paper. The superim- 
posed false-color pattern represents the value of the Ez field 
in the {x = 0, y, z) plane (the blue and red colors respectively 
correspond to the maximum and minimum values of the field) . 



We consider a cylindrical-symmetric pillbox cavity op- 
erating on the TMiio at the zero-crossing phase. The 
corresponding electromagnetic field components in the 
paraxial approximation are given by 

Ez{x,y, z,t) — EQXCOs{iut) 



By{x,y,z,t) 



E'„ 



^ sin(a;i). 



(29) 
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where E'q = dE^/dx ~ {2tt/X)Eq (where Eq is the peak 
electric field) and cj = 27r/ (/ is the operating frequency). 
The thick-lens transfer matrix in {x,x' ,z,5) is given by 
Eq. [2] where k = ^ ^^ . The cavity used in our simu- 
lations is shown in Fig. [1^1 It operates at / = 1.3 GHz 
and consists of three cells of length (1/2 — 1 — 1/2) x A/2. 
Two end pipes were added in the electromagnetic model 
to properly account for the fringe fields of the cavity. 

Considering the cavity to consists of three uncoupled 
pillbox resonators of length A/4, A/2 and A/4 we compute 
the transfer matrix to be 



R 



Ru Rl2 Rl5 ^16 

R21 R22 R25 R2& 

R5I -R52 ^55 -R56 

. Rqi i?62 ^65 -^66 



1 ic 

1 



K ^ 



2 
K 

1 

23k^A 
128 



(30) 



where Lc is the cavity length including the outer pipes; 
see Fig. [HI 

To numerically evaluate the first order transfer matrix 
of the cavity, we track seven macroparticles in the cavity 
showed in Fig. [TB] Six of the macroparticles are displaced 
by an infinitesimal small amount (5X* with respect to 
the reference macroparticle along one of the phase space 
coordinates. This initial displacement is mapped down- 
stream of the cavity as a displacement (5X'"' with respect 
to the final position of the reference macroparticle. The 
initial and final displacement are related via 



vides the transfer matrix elements. 
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FIG. 17: Comparison of the analytical (solid lines) and nu- 
merical (squares) calculations of the /i-dependent elements of 
transfer matrix associated to the deflecting cavity considered 
in this paper. 



The K-dcpcndcnt elements of the numerically- 
computed transfer matrix are compared with those ana- 
lytically calculated using Eq. [30] in Fig. [171 The agree- 
ment between the numerical and analytical models is 
excellent; the relative discrepancy is < 5% over k, S 
Inversion of the latter system of equations directly pro- [~5, 5] m^^. 



SX! -j = y^RijS^i 



(31) 
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